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1. Introduction 

Lattice simulations performed by different groups involve different choices both at the level 
of formalism (lattice actions, number of sea flavours etc.) and at the level of resources (lattice 
volumes, quark masses etc.). Often this amounts to making different compromises which in turn 
introduce different systematic effects; thus not all lattice results of a given quantity are directly 
comparable. The Flavour Lattice Averaging Group (FLAG) aims to answer, in a way which is 
readily accessible to non-experts, the question: When it comes to lattice quantiities of relevance to 
Flavour Physics, what is currently their “best lattice value”? Two editions of the FLAG review have 
appeared so far [1, 2], to which we will be referring as FLAG-1 and FLAG-2 for brevity. Currently 
the FLAG-3 report is under preparation ^ and should appear in early 2016. 

In the present summary we will only show a subset of the results reviewed in FLAG-2 [2], 
selecting those which are closely related to the topic of this conference (Flavour Physics & CP- 
violation). The FLAG-2 compilation was based on results published in refereed journals by the 
30th of November 2013, so they are somewhat outdated, while those of FLAG-3 is still in the 
works. As a compromise we show FLAG-2 results for //) 5 , the D- and 5-meson semileptonic 
decay form factors, and the corresponding CKM matrix elements, since their update is still in 
a rather preparatory phase. We present preliminary FLAG-3 results for (and also briefly for 
fo), the ^-meson semileptonic decay form factor, the corresponding CKM matrix elements and Bk, 
since their update is in an advanced stage. We also review two topics which go beyond the scope 
of FLAG: (i) an interesting analysis [3] of the Standard Model (SM) prediction of £k, based on the 
FLAG-2 estimate of Bk\ (h) the kaon oscillation 5-parameters beyond the SM, which appeared in 
Refs. [4, 5]. 

FLAG analyses are based on a critical review of lattice results. A number of criteria have been 
fixed, providing compact information on the quality of a computation. These are typically related to 
the quality of continuum and chiral extrapolations, finite volume effects, etc. The FLAG-2 rating is 
organized in Tables, in which different results are colour-coded as follows: ★ means that the rele¬ 
vant systematic error has been estimated in a satisfactory manner and is under control; o means that 
a reasonable attempt at estimating the systematic error has been made, which however can be im¬ 
proved; ■ means that there was no attempt or that the attempt at controlling a systematic error was 
unsatisfactory. In the latter case the result is dropped from FLAG averages or estimates. The Tables 
are not presented here for lack of space. FLAG-3 adopted a similar rating, with somewhat differ¬ 
ently phrased criteria. Several other issues are addressed in Ref. [2]; e.g. how to average different 
results; how to arrive at a final estimate if an average is not possible; how to combine/correlate 
errors; how (not) to take conference proceedings into account. The interested reader should consult 
the introductory sections of Ref. [2] for detailed explanations. For our purposes it suffices to keep 
in mind the following two FLAG policies: 

Firstly, we recall that lattice simulations are carried out with a fixed number Nf of dynamical 
(sea) quarks. So far results exist for A/ = 0 (quenched), Nf — 2 (two degenerate light quarks) Nf = 
2 -h 1 (including a heavier strange quark) and Nf — 2+\ + \ (as before, plus charm). Quenched 
results are outdated and are thus omitted, with the exception of those related to see Fig. 1. 

^For the composition of the collaboration, past and present, see the FLAG website: 
http://itpwiki.unibe.ch/flag 
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FLAG averages or estimates are quoted separately for each Nf \ results referring to different Nf are 
never averaged, as they correspond to different field theoretic approximations of QCD. 

Secondly, our Figures are organized as follows: results are clearly separated according to their 
Nf value. Each datapoint is accompanied by an acronym, indicating the relevant reference, cited in 
Ref. [2]. Lack of space prevents us from citing all these works here; however we do cite, whenever 
possible, the papers containing the results on which FLAG averages/estimates are based. The plots 
are colour coded as follows: 

■ indicates a FLAG average or estimate; these results are also highlighted by a gray vertical band; 

■ refers to data that have passed our quality criteria and are used in the determination of the FLAG 
results (i.e. the points indicated by ■); 

□ indicates results with good control of the systematics (i.e. no red tags in the Tables of Ref. [2], 
which are left out of the average for some reason; e.g. not published in peer-reviewed journals, 
superseded by later results of the same collaboration, or simulations suffering from some uncon¬ 
trolled effect which is not colour coded in the Tables)^ 

□ refers to results that are not included in the average because they do not pass the rating criteria 
(i.e. they have at least one ■ in the Tables or do not pass for some other shortcoming); 

• indicates non lattice results, presented for comparison. 

A good example of the effects of our rating in our final result is provided by the analysis of as 
in Ref. [2]. In Fig. 1 we show the strong coupling for 5 flavours, at the scale Mz in the MS scheme. 
Preferring to err on the side of caution, the majority of FLAG-2 members opted for a conservative 
error estimate, based on Refs. [6, 7, 8, 9], citing: 

ag(Mz) =0.1184(12). (1.1) 


This is fully compatible with the PDG average, obtained by excluding lattice results: o^(Mz) = 
0.1183(12). Note that the error of the FLAG lattice estimate is as big as that of the above PDG 
non-lattice result. 


2. Light Flavour Physics 

We summarize state of the art lattice calculations of the leptonic kaon and pion decay constants 
and the kaon semileptonic decay form factor and provide an analysis in view of the Standard Model. 
The pion decay constant is defined by the matrix element of the axial current = djfxYsu between 
vacuum and pion states: 

(0|A^|7r+(p)) = ipnfj,+ (2.1) 

In this normalization, fjj^± 130 MeV. Analogous expressions hold for fK+ and (in the following 
sections) the decay constants of charmed and bottom mesons. The kaon semileptonic decay form 
factor /+(^^) represents one of the form factors relevant for the semileptonic decay -A 
which depends on the momentum transfer (squared) between the two mesons. What matters 
here is its value at = 0. 

^In FLAG-3 this square is filled with ligh green; cf. Figs. 2,3,7. 
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Figure 1: the coupling constant in the MS scheme at the Z mass. For the reference labels see 

Ref. [2]. The results labeled A/ = 0,2 use estimates for Nf = 3, obtained by first extrapolating in Nf from 
A/ = 0,2 results. Since this is not a theoretically justified procedure, these are not included in our final 
estimate and are thus given a red symbol. The black triangle indicates the outcome of the PDG analysis 
excluding lattice results. 


All lattice results for the form factor /+(0) and many available results for the ratio of decay 
constants have been computed in isospin-symmetric QCD. The reason for this unphysical param¬ 
eter choice is that there are only few lattice simulations of SU(2) isospin-breaking effects. When 
combining lattice data with experimental results, we take into account the strong SU(2) isospin 
correction, either obtained in lattice calculations or estimated by using chiral perturbation theory, 
both for the kaon leptonic decay constant and for the ratio Here and fjj;± are the 

isospin-broken decay constants in QCD (the electromagnetic effects have already been subtracted 
in the experimental analysis using chiral perturbation theory). The plots in Figure 2 illustrate the 
preliminary FLAG-3 compilation of data for /+(0) and fj^±/ fj^±. The lattice data for the latter 
quantity are largely consistent even when comparing simulations with different A/, while in the 
case of /+ (0) a slight tendency to get higher values for increasing A/ seems to be visible, not ex¬ 
ceeding one standard deviation. The A/ = 2 + 1 FLAG-average for /+ (0) is based on FNAL/MILC 
12 [10] and RBC/UKQCD 15 [11], which we consider uncorrelated, while for A/ = 2 + 1 + 1 and 
A/ = 2 we consider directly the FNAL/MILC 14 [12] and ETM 09A [13] results, respectively: 


/+(0) = 0.9704(24) (22), 

(iV/ = 2+l + l) 

(2.2) 

/+(0) = 0.9677(37), 

(Nf = 2+1) 

(2.3) 

/+(0) = 0.9560(57) (62). 

of 

11 

(2.4) 


The brackets in the first and third lines indicate the statistical and systematic errors respectively. 
We stress that the results (2.2) and (2.3) include simulations with physical light quark masses. Note 
that the Nf = 2 result remains unchanged from FLAG-2. 
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Figure 2: Comparison of lattice results (squares) for /+ (0) and /^± / f^± with various model estimates based 
on ;^PT (blue circles). The ratio f^± /is obtained in pure QCD including the SU(2) isospin breaking 
correction. For the reference labels see Ref. [2]. 


For the deacy constant ratio j we quote: 


/^±//;r± = 1-193(3), {Nf = 2+\ + \) 

fK^/U = 1.192(5), (Nf = 2+1) 

fK±/A± = l.205{6){n), {Nf = 2) 


(2.5) 

( 2 . 6 ) 
(2.7) 


where in the last row the first error is statistical and the second systematic. The first row is a 
preliminary update of FLAG-2, whereas the last two have remained unchanged. 

Precision experimental data on kaon decays very accurately determine the product \Vus\f+{0) 
and the ratio |K./14j|/jf±//;r± [14]: 


|l4,|/+(0) = 0.2163(5) 



Iei 

A=*= 


= 0.2758(5) 


( 2 . 8 ) 


Combining the above with our lattice estimates for /+(0) and we obtain the angles \Vud\ 

and \Vus\, thus testing first row unitarity of the CKM matrix: 


\Vu\^=\Vud?+\yus?+\Vub?=l 


(2.9) 


The tiny contribution from \Vub\ is known much better than needed in the present context: \Vub\ = 
4.15(49) • 10“^ [15]. Evidence for the validity of the relation (2.9) is shown in Figure 3 (left). 
From our results for /+(0) and fK^/fn^ we obtain a corresponding range for the CKM matrix 
elements \Vud\ and \Vus\, using the relations (2.8). Consider first the results for A/ = 2 + 1 + 1. 
The range for /+(0) is mapped into the interval \Vus\ = 0.2229(9), depicted as a horizontal red 
band in Figure 3 (left), while the one for is converted into =0.2311(7), shown 

as a tilted red band. The red ellipse is the intersection of these two bands and represents the 
68% likelihood contour, obtained by treating the above two results as independent measurements. 
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Repeating the exercise for A// = 2 + 1 and Nf = 2 leads to the green and blue ellipses, respectively. 
The correlation between \Vud \ and \Vus \ imposed by the unitarity of the CKM matrix is indicated by 
a dotted arc (in view of the uncertainty in \Vub\, this is really a band of finite width, but the effect is 
too small to be seen here). The plot shows that there is a slight tension with unitarity in the data for 
A// = 2 + 1 + 1: Numerically, the outcome for the sum of the squares of the first row of the CKM 
matrix reads |14p = 0.980(10), which deviates from unity at the level of two standard deviations. 
In spite of this tension, it is fair to say that at this level the Standard Model passes a nontrivial test 
that exclusively involves lattice data and well-established kaon decay branching ratios. If we use 
the P decay value of \Vud\ = 0.97425(22) quoted in Ref. [16], the test sharpens considerably: the 
lattice result (2.2)) for /+(0) leads to |I4p = 0.9989(8), while the one for fx^/ (2.5) implies 
|I4p = 0.9999(7), thus confirming CKM unitarity at the permille level. The situation is similar for 
A/ = 2 + 1: |Kp = 0.983(12) from the lattice data alone. Using the j8 decay value of \Vud\ again 
the test sharpens, giving | 14 p = 0.9991(9) from the lattice value for /+(0) and | 14 p = 1.0000(6) 
from ff^± // 7 r±- Repeating the analysis for Nf — 2, we find |14p = 1.029(35) with the lattice data 
alone, which is fully compatible with unity. When using the nuclear /3 decay value of \Vud\, we 
obtain |14p = 1.0004(10) from/+(0) and |14p = 0.9989(16) /fj^±. 



Vud 


FLAG2015 



Figure 3: Left: The plot compares the information for |14j|, 114^1 obtained on the lattice with the exper¬ 
imental result extracted from nuclear j3 transitions. The dotted arc indicates the correlation between \Vud\ 
and 114^1 that follows if the three-flavour CKM-matrix is unitary. Right: Results for 114^1 and \Vud\ that 
follow from the lattice data for /+(0) (triangles) and fx±/fji± (squares), on the basis of the assumption 
that the CKM matrix is unitary. The black squares and the grey bands represent our estimates, obtained by 
combining these two different ways of measuring and \ Vud \ on a lattice. For comparison, the figure also 
indicates the results obtained if the data on nuclear j8 decay and T decay are analyzed within the Standard 
Model. For the reference labels see Ref. [2]. 

The Standard Model implies that the CKM matrix is unitary. The precise experimental con¬ 
straints quoted in (2.8) and the unitarity condition (2.9) then reduce the four quantities \Vud\, \Vus\, 
/+(0), to a single unknown: any one of these determines the other three within narrow 
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uncertainties. Fig. 3 (right) shows that the results obtained for \Vus\ and \Vud\ from the data on 
fK^/fn^ (squares) are quite consistent with the determinations via /+(0) (triangles). In order to 
calculate the corresponding average values, we rely on those results that have passed the various 
rating criteria; for details see Ref. [2] and the forthcoming FLAG-3 update. The comparison shows 
that the lattice result for \Vud \ not only agrees very well with the totally independent determination 
based on nuclear j8 transitions, but is also remarkably precise. On the other hand, the values of 
\Vud\ which follow from the T decay data if the Standard Model is assumed to be valid, are not in 
good agreement with the lattice results. The disagreement is reduced considerably if the analysis 
is supplemented with experimental results on electroproduction. The preliminary FLAG-3 results 
for \Vus\ and \Vud\ (assuming unitarity of the CKM matrix) are 


\Vus\ 

= 0.2248(9), 

iKrfI =0.97440(21), 

(V/ = 2 + l + l) 

(2.10) 

IK.1 

= 0.2245(10), 

\Vud\ =0.97447(23), 

(V/ = 2+l) 

(2.11) 

|v«.| 

= 0.2253(21), 

\Vud\ =0.97427(49). 

iNf = 2) 

(2.12) 


We close this section by quoting, without details, results for the decay contants and f^. 
This is possible only for works which have not used fji for “setting the scale” (i.e. for the deter¬ 
mination of the finite UV cutoff in the simulations). For the pion decay constant we have only the 
Flag-2 [2] estimate for A// = 2-h 1; = 130.2(1.4) MeV, to be compared with the PDG [15] value 

= 130.41(20) MeV. For we quote the FLAG-3 preliminary values: 

fK± = 155.6 (0.4) MeV (V/ = 2+l + l), (2.13) 

fK± = 155.9 (0.9) MeV (V/ = 2-h 1), 

= 157.5 (2.4) MeV (A^/= 2), 

to be compared with the PDG [15] result = 156.2(7) MeV. 

3. Charm Physics 

Leptonic and semileptonic decays of charmed D- and mesons are sensitive probes of c^d 
and c ^ s quark flavour-changing transitions. They enable the determination of the CKM matrix 
elements 114^1 and \Vcs \ (within the Standard Model) and a precise test of the unitarity of the second 
row of the CKM matrix. Here we summarize the status of lattice-QCD calculations of the charmed 
leptonic decay constants and semileptonic form factors ^(0), as reported in FLAG-2 [2] and 
shown in Figure 4. Our final averages are: 

//) = 209.2(3.3) MeV,//)^ = 248.6(2.7) MeV, ^ = 1.187(12) (A^y = 2+1) (3.1) 

JD 

/o = 208(7) MeV,/o, = 250(7) MeV, ^ = 1.20(2) iNf = 2) (3.2) 

Id 

The above Nf = 2 is the result of the only reference in the literature at the time, ETM 13B [17]. The 
Nf = 2+1 estimates for fo and the SU(3) breaking ratio foJfD come from HPCQD 12A [18] and 
FNAL/MILC 11 [ 1 9] , while the result for fo, comes from HPQCD 1OA [20] and FNAL/MILC 11 [ 1 9] . 
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Figure 4: Decay constants of the D and Ds mesons. Errors in FNAL/MILC 13 are smaller than the symbols. 
For the reference labels see Ref. [2]. 


Because FNAL/MILC and HPQCD use a largely overlapping set of configurations, we treat the sta¬ 
tistical errors as 100% correlated. A preliminary FLAG-3 analysis gives fo, = 249.8(2.3), upon 
including ;tQCD 14sea [21]. At the time of writing the FLAG-2 report the first A/ = 2 + 1 + 1 
computations had only appeared in conference proceedings, so no averages had been worked out. 
Based on Refs. [22] (FNAL/MILC 14wgs) and [23] (ETM 14poa), a preliminary analysis gives 
fD = 212.15(1.12) MeV, fo, = 248.83(1.27) MeV md/oj/o = 1.1716(32). We stress that since 
the accuracy of the lattice determinations of the D meson decay constant is rapidly improving, 
it will become important in the future, especially when comparing to experimental numbers, to 
distinguish between //)+ and the average of //)+ and /^o. 

We now turn to the form factors for semileptonic D Ttiv and D Kiv decays. In practice, 
most lattice-QCD calculations focus on providing the vector form factor at a single value of the 
momentum transfer, /+(^^ = 0), which is sufficient to obtain \Vcd\ and \Vcs\- Because the decay 
rate cannot be measured directly at zero momentum transfer, comparison of these lattice-QCD 
results with experiment requires a slight extrapolation of the experimental measurement. FLAG-2 
did not quote any averages/estimates for A/ = 2 and A/ = 2 + 1 + 1, as at the time they were at a 
preliminary stage. For A/ = 2 +1 there was only one simulation satisfying all of our quality criteria 
(HPQCD lOB [24]), HPQCD 11 [25]) from which we quote: 

/^^(o) =0.666(29), /^^(O) =0.747(19) (iV/ = 2+l). (3.3) 

The results for the Di^s) meson decay constants (related to the branching ratio for leptonic de¬ 
cays), as well as the semileptonic decays’ form factors, lead to determinations of the CKM matrix 
elements \Vcd\ and m the Standard Model. For the leptonic decays, we use the latest experi¬ 
mental averages from Rosner and Stone for the Particle Data Group [26] (where electromagnetic 
corrections of ^ 1% have been removed): 

=46.40(1.98) MeV, = 253.1(5.3) MeV. (3.4) 
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Figure 5: Left: Comparison of determinations of \Vcd\ and \Vcs\ obtained from lattice methods with non¬ 
lattice determinations and the Standard Model prediction based on CKM unitarity. When two references are 
listed on a single row, the first corresponds to the lattice input for \Vcd\ and the second to that for IVc^j. The 
results denoted by squares are from leptonic decays, while those denoted by triangles are from semileptonic 
decays. For the reference labels see Ref. [2]. Right: Comparison of determinations of \Vub\ obtained 
from lattice methods with non-lattice determinations based on inclusive semileptonic B decays. The results 
denoted by squares are from leptonic decays, while those denoted by triangles are from semileptonic decays. 

We combine these with the average values of fo and in Eqs. (3.2) and (3.1), to obtain \Vcd\ and 
\Vcs\- For the semileptonic decays, we use the latest experimental averages from the Heavy Flavour 
Averaging Group (HFAG) [33]. 

/?y0)|Krf| = 0.146(3), = 0.728(5). (3.5) 

and /+^(0) and /+^(0), of Eq. (3.3), to obtain our preferred values for \Vcd\ and \Vcs\- Fig. 5 dis¬ 
plays the results for 114^1 and \Vcs\ from leptonic and semileptonic decays, and compares them to 
determinations from neutrino scattering (for 114^1 only) [15] and CKM unitarity [26]. The determi¬ 
nations of \Vcd \ agree within uncertainties. The determination of \Vcs\ from Nf = 2+1 lattice-QCD 
calculations of leptonic decays is noticeably larger than that from both semileptonic decays and 
CKM unitarity. The disagreement between \Vcs\ from leptonic and semileptonic decays is slight 
(only 1.2(7 assunfing no correlations), but the disagreement between \Vcs \ from leptonic decays and 
CKM unitarity is larger at 1.9(7. This tension is driven primarily by the HPQCD lOA [20] calcula¬ 
tion of but we note that the Nf = 2 calculation of (ETM13B [17]) leads to the same high 
central value of | 14 .yUttst with larger uncertainties. 

The A/ = 2 + 1 averages for \Vcd\ and \Vcs\ in Fig. 5 take correlations of different results into 
account. Omitting details which the reader may find in Ref. [2] , we simply state that whenever 
there are correlations, they are assumed to be 100%. We obtain 

= 0.2191(83), |K,| = 0.996(21), (Nf = 2+1) (3.6) 

where the errors include both theoretical and experimental uncertainties, and the error on \Vcs\ has 
been increased by = 1.03. 
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Using the determinations of \Vcd\ and \Vcs \ in Eq. (3.6), we can test the unitarity of the second 
row of the CKM matrix. We obtain 

\Vcdf + - 1 = 0.04(6) (3.7) 

which agrees with the Standard Model at the percent level. Given the current level of precision, 
this result does not depend on the value used for \Vcb\, which is of 

4. Bottom Physics 

Lattice-QCD calculations of Z^-quarks have an added complication not present for charm and 
light quarks: at the lattice spacings that are currently used in numerical simulations, the Z^-quark 
mass is of order one in lattice units. As a direct treatment of Z^-quarks with the fermion actions 
commonly used for light quarks would result in large cutoff effects, all current lattice-QCD calcu¬ 
lations of Z 7 -quark quantities make use of effective field theory at some stage. The two most widely 
used general approaches are (i) direct application of effective field theory treatments such as HQET 
or NRQCD, which allow for a systematic expansion in 1/m^; or (ii) the interpretation of a relativis¬ 
tic quark action in a manner suitable for heavy quarks using an extended Symanzik improvement 
program to suppress cutoff errors. This introduces new problems (matching of HQET to QCD, 
renormalization, control of discretization effects...). 

In the present summary we will only quote results of the leptonic decay constants fs and 
and the CKM matrix element \Vub\, once more leaving details and other results (on semileptonic 
form factors of heavy-to-light and heavy-to-charm decays, as well as neutral 5-meson mixing ma¬ 
trix elements) to Ref. [2]. The leptonic decay constants, shown in Eig. 6 , lead to the following 
ELAG-2 estimates: 

fs = 189(8) MeV, = 228(8) MeV, ^ = 1.206(24) {Nf = 2) (4.1) 

JB 

fB = 190.5(4.2) MeV, fe, = 227.7(4.5) MeV, ^ = 1.202(22) (V/ = 2+1) (4.2) 

JB 

/g = 186(4) MeV,/ b^= 224(5) MeV, ^ = 1.205(7) (V/ = 2+l + l) (4.3) 

JB 

Most results, obtained with degenerate light quarks, refer to average decay constants for 5+ and 

Some collaborations (ENAL/MILC, HPQCD) have started giving distinct results (they differ 
by about 2%). As errors decrease with time, collaborations should start giving 5+ and results 
separately. The preliminary ELAG-3 update shows no significant changes for Nf = 2 and Nf = 
2+1, while the Nf = 2 + \ + \ numbers remain unaltered. 

We now use the lattice-determined Standard Model transition amplitudes for leptonic and 
semileptonic 5-meson decays to obtain exclusive determinations of the CKM matrix element |14^|. 
The branching fraction for the decay 5 ^ has been measured by the Babar [27, 28] and 
Belle [29, 30] experiments. Combining the averaged data by the two experiments and the lattice 
ELAG2 averages, ref. [2] quotes: 


|V„fe|=4.21(53)(18)x 10-3 

V/ = 2, 

(4.4) 

|K;,|=4.18(52)(9)x10-3 

iV/ = 2+l, 

(4.5) 

|Kfc|=4.28(53)(9)x 10-3 

IV/-= 2+1 + 1. 

(4.6) 
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Figure 6: Decay constants of the B and Bs mesons. For the reference labels see Ref. [2]. 

In semileptonic B^ decays we use as experimental input the Belle [31] and Babar [32] 

resuts, combined with lattice form factor estimates for A// = 2 + 1, quoting: 

global lattice + Babar: \Vub\ = 3.37(21) x IQ-^, Nf = 2+\, (4.7) 

global lattice + Belle: \Vub\ = 3.47(22) x IQ-^, = 2 + 1. (4.8) 

We do not quote a result for a combined lattice + Babar + Belle fit, since we are unable to properly 
take into account possible correlations between experimental results. 

Our results for \Vub\ are summarized in Figure 5 (right), where we also show the inclusive de¬ 
terminations from HFAG [33] for comparison. The spread of values for \Vub\ does not yield a clear 
picture. We observe the well-known ^3(7 tension between determinations of \Vub\ from exclu¬ 
sive (lattice) and inclusive (HFAG) semileptonic decays. The determination of \Vub\ from leptonic 
B ^ TV decay lies between the inclusive and exclusive determinations, but the experimental errors 
in BR{B tv) are so large that it agrees with both within ^ 1.5(7. The exclusive determination of 
\Vub\ will improve in the next few years with better lattice-QCD calculations of the B ^ niv form 
factor, while the improvement in \Vub\ from B ^ TV decays will have to wait longer for the Belle II 
experiment, which aims to begin running in 2016. 

5. Neutral AT-meson oscillations and Bk 

Indirect CP-violation in Kaons is measured by the parameter Ek. Followig Ref. [3] we express 
it here as 

Ek — Vlsinoi^^BKXsD + <^o + <^ld^ H- (5.1) 

where some tiny correction terms have been dropped. This expression is valid in QCD with three 
light flavours; i.e. an effective theory, in which charm and heavier degrees of freedom have been 
integrated out. The factor Cg = [G^Fjm^oM^]/[6v^;r^AM^] is a known quantity, and Bk is the 
bag parameter, which incorporates low-energy QCD effects. It is defined through the hadronic 
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matrix element of a four-fermion, dimension-six operator normalized by factors of the Kaon 
mass and decay constant: 




{K^\Qf=\^)\K^) 

^fK^K 


(5.2) 


Note that 5 k(m) is a renormalized quantity (thus the subscript in which depends on a chosen 
renormalization scheme and scale ji. Its multiplication by a renormalization group evolution func¬ 
tion, known in NLO perturbation theory, results in the scale independent (renormalization group 
invariant) quantity Bk which appears in eq. (5.1). 5 k(m) subsequently Bk) is now known from 
several lattice computations, which have been reviewed by FLAG. The short-distance contribution 
is 


XsD = ri^^\Vcb\^ 


~p)%5o(x^)(l + r) -h 



{ Bct^O (-^c 5 ) Bcc^O (-^c ) } 


(5.3) 


where So{xc^t) and 5o(xc,x^) are the Inami-Limfunctions (with Xc^t = ^ct/^w^ Vtt^Vct^Vcc are 
factors known perturbatively to NLO, NNLO and NNLO respectively; the quantity r is a known 
function of 5o’s and T]’s. Long-distance effects from the absorptive part are estimated to have 
a —7% contribution [34], while long-distance effects <^ld from the dispersive part contribute only 
2% [35] and will be neglected in the following discussion. For more details see Ref. [3] . 

In Fig. 7 we present the preliminary compilation of results of FLAG-3. For A// = 2 -h 1, the 
result is^ 

Bk = 0.7627(97) Nf = 2+1 (5.4) 

This is compatible with the = 2 -h 1 value Bk = 0.7661 (99) quoted in FLAG-2 [2]. 

The FLAG-2 result has been used in [3] for an interesting phenomenological analysis: Eqs. (5.1) 
and (5.3) are used for a theoretical (Standard Model) prediction of £k and the value is com¬ 
pared to the known experimental result. For the Wolfenstein parameters p and fj they prefer 
the so-called Angle-Only-Fit results of Ref. [36] to the more standard UTfit / CKMfitter ones, 
because the latter contain unwanted dependence on Bk, \Vcb\ and £k. They use \Vus\ ~ A from 
and K /3 decays and 114^1 AX^. From Ref. [37] on inclusive B XJv and B XsY de¬ 
cays they take 114^1 = 42.21(78) x 10“^. From Ref. [38] on exclusive B D^lv decays they use 
|K:Z7H39.04(49)(53)(19) x IQ-^ These result to \£^\ = 1.58(18) x 10“^ (exclusive) and \£^\ = 
2.13(23) X 10“^ (inclusive). Comparison with the experimental value \£k^\ = 2.228(11) x 10“^ 
shows a stress of A8 k = I^k^I — = 3.6(2)<7 in the exclusive case, which is absent in the 

inclusive one (A8 k = 0.44(24)o'). The former case requires further investigation. 

A second interesting observation made in [3] concerns the error budget they report for \£^\, 
obtained with exclusive \Vcb\ and the FLAG-2 value of Bk- Its error is dominated by |t4^| (nearly 
41%), which is not surprising, given that \Vcb\ enters with a fourth power in Eq. (5.3). In compari¬ 
son, the contribution of the Bk error is a mere 1.6%, which indicates that lattice results nowadays 
are determined with errors which are clearly subdominant in certain phenomenological analyses. 

We close this section with a brief discussion of possible New Physics (NP) effects related to 
neutral Kaon oscillations, as recently analyzed in a model-independent way. A generalization of 

^For comparison we also give Bk = 0.727(25) for Nf = 2 and Bk = 0.717(24) for A/y = 2 -h 1 -h 1. 
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Figure 7: Left: Lattice results for the renormalisation group invariant ^^^-parameter. For the reference 
labels see Ref. [2]. Right (top): The lower bounds on the New Physics (NP) scale A, for generic NP flavor 
structure, obtained from the Nf = 2 simulation of Ref. [4] (brown bars), compared to the bounds obtained 
form the quenched analysis of Ref. [39] yellow bars). Right (bottom): A compilation of meson bag- 
parameters / = 2,..., 5. From top to bottom data have been taken from Refs. [4, 40, 41, 5]. 


the effective A5 = 2 Hamiltonian is assumed, which contains operators beyond the one arising in 
the SM ; the amplitude is: 


5 

= Ci(^Vi|^°) + • (5-5) 

i^2 

The first term on the rhs is the SM contribution"^ and 02,.. .,5 are four-fermion operators of dimension- 
6, which parametrize NP effects. In analogy to Bk = Bi, one defines 52, -,5- The Wilson coedffi- 
cients C 2,...,5 are parametrized as Ci = [FiLi]/A^, where Fi stads for NP couplings, Li are coupling 
dependent loop factors and A is the scale of NP. Assuming, for a generic strongly interacting theory 
with an unconstrained flavor structure, F; ^ L/ ^ ^(1), a generalized UT-fit analysis produces lower 
bounds for A; these depend very strongly (by several orders of magnitude!) on this assumption. To 
avoid accidental cancellations each contribution is analyzed separately. Fig. 7 (right, top) compares 
these lower bounds on A resulting from two lattice simulations. They differ by the number of dy¬ 
namical quark flavours in the sea, one being quenched {Nf — 0) and the other having Nf = 2. More 
importantly, they also differ by the accuracy with which the Bi parameters have been determined 
(20% — 23% for Nf — Q and 3% — 6% for Nf = 2). The SM bound turns out to be several orders of 
magnitude weaker than those arising from NP operators. 

"^i.e. Oi = and {K^\0\ is the numerator of Bk in Eq. (5.2). 
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This leads to the need of a more stringent determination of the NP 5-parameters. Several 
groups are producing results, which will be analyzed in FLAG-3. As a foretaste we show the 
recent comparison by Ref. [5] in Fig. 7 (right, bottom). The situation is still unsettled, especially 
as far as 54 and are concerned, in view of the different souces of systematic error characterizing 
the data. 
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